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Abstract 

Rarefied  gas  flows  through  two-dimensional  fi¬ 
nite  length  capillaries  are  studied  numerically 
with  the  goal  of  performance  optimization  of  a 
nanomembrane-based  Knudsen  compressor.  Both 
density  and  temperature  gradient  driven  flows  are 
considered  in  a  wide  range  of  Knudsen  numbers 
from  0.05  to  75.  The  length-to-height  ratios  from 
5  to  30  were  examined.  Three  different  wall  tem¬ 
perature  distributions  were  considered,  namely,  lin¬ 
ear,  step-wise,  and  a  non-monotonic  profile  typical 
for  a  radiantly  heated  Knudsen  compressor  mem¬ 
brane.  The  direct  simulation  Monte  Carlo  method 
was  used  for  shorter  channels,  and  a  discrete  or¬ 
dinate  method  for  the  solution  of  the  BGK  model 
equation  is  utilized  for  a  longer  channel.  The  mech¬ 
anism  for  the  Knudsen  minimum  in  flow  conduc¬ 
tance  is  quantitatively  explained,  the  short  channel 
end  effects  are  characterized,  and  the  sensitivity  of 
the  mass  flow  rate  to  a  non-monotonic  temperature 
distribution  is  shown. 

1  Introduction 

Lab  on  a  chip  technology  is  drawing  the  attention  of 
scientists  from  many  disciplines.  Recent  advances 
in  MEMS  manufacturing  have  made  it  possible  to 
construct  microscale  analytical  sensors  such  as  in¬ 
tegrated  gas  chromatography  systems,  miniature 
spectrometers,  and  mass  spectrometers.  The  minia¬ 
turized  detection  devices  will  require  microscale 
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roughing  pumps  to  provide  sensor  elements  with  gas 
samples  at  necessary  environmental  conditions.  A 
pumping  mechanism  that  can  be  exploited  at  mi¬ 
croscale  is  thermal  transpiration,  a  rarefied  gas  ef¬ 
fect  that  drives  the  gas  flow  along  the  temperature 
gradient  in  a  tube  or  channel.  The  main  goal  of 
this  paper  is  the  numerical  study  of  thermal  tran¬ 
spiration  flows  in  short  microchannels  to  aid  in  the 
performance  optimization  of  a  transpiration  based 
microscale  roughing  pump  known  as  the  Knudsen 
compressor  [1], 

In  early  1900s  M.  Knudsen  built  and  studied  the 
first  transpiration  based  compressor,  consisting  of  a 
series  of  differentially  heated  and  cooled  capillaries. 
Each  stage  of  the  compressor  had  a  capillary  sec¬ 
tion  where  the  wall  temperature  increases,  causing 
a  thermomolecular  pressure  build-up  at  the  high- 
temperature  end  of  the  capillary.  The  capillaries 
are  followed  by  a  connector  stage  of  a  significantly 
larger  cross-sectional  area,  where  the  pressure  is 
almost  constant  and  temperature  is  decreased  to 
its  original  value  at  the  beginning  of  the  stage. 
The  modern  version  of  Knudsen  compressor  was 
suggested  by  Pham-Van-Diep  et  al[ 2]  and  demon¬ 
strated  by  Vargo  and  Muntz  [3].  The  most  critical 
element  of  the  Knudsen  compressor  developed  at 
USC  is  the  thermal  transpiration  membrane  made 
of  porous  materials,  such  as  aerogel,  with  pore  di¬ 
ameters  on  the  order  of  the  mean  free  path  of  the 
gas  at  standard  atmospheric  conditions.  The  tem¬ 
perature  gradient  is  maintained  across  the  transpi¬ 
ration  membrane  by  resistive  or  radiant  heating. 
Recently,  a  single-chip  micromachined  implementa¬ 
tion  of  a  Knudsen  pump  was  reported  in  Ref.  [4], 
which  can  evacuate  a  cavity  to  0.46  atm  while  op¬ 
erating  at  atmospheric  pressure  and  using  80mW 
input  power. 

Knudsen  compressor  has  several  attractive  fea- 
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tures,  such  as  no  moving  parts  and  no  fluids  or  lubri¬ 
cants,  and  therefore,  can  have  a  number  of  promis¬ 
ing  applications  provided  that  its  performance  can 
be  optimized  to  meet  the  existing  pumping  re¬ 
quirements.  The  two  most  important  performance 
parameters  of  the  compressor  are  the  energy  use 
and  volume  per  unit  throughput.  To  accomplish 
minimization  of  these  performance  parameters  the 
Knudsen  compressor  must  also  be  designed  to  max¬ 
imize  the  gas  conductance  through  the  device  for  a 
given  temperature  change.  The  required  flow  rate 
in  the  Knudsen  compressor  can  be  obtained  by  ad¬ 
justing  the  cross-sectional  area  of  the  stages  and 
the  temperature  distribution  along  the  transpira¬ 
tion  membrane  at  different  stages.  The  main  el¬ 
ement  of  the  transpiration  based  pump  is  rarefied 
gas  flow  through  a  capillary. 

Rarefied  gas  flows  of  different  gases  through  long 
tubes  and  channels  were  a  topic  of  an  extensive 
experimental  and  theoretical  research  over  several 
decades,  and  a  comprehensive  review  of  the  nu¬ 
merical  data  and  analytical  results  for  rarefied  gas 
flows  in  capillaries  was  performed  by  Sharipov  and 
Seleznev[5]  (see  also  references  therein).  There  is 
however  a  lack  of  reference  data  for  nonisothermal 
flow  through  a  capillary  of  a  finite  length,  especially 
for  arbitrary  pressure  and  temperature  drops. 

The  present  numerical  modeling  of  thermal  tran¬ 
spiration  flows  in  microchannels  addresses  sev¬ 
eral  questions  important  for  design  optimization  of 
Knudsen  compressors,  namely:  the  influence  of  dif¬ 
ferent  temperature  internal  distributions  across  the 
membrane  on  the  flow  field  inside  the  capillaries  and 
the  total  transpiration  mass  flow  rates,  the  impact 
of  the  rarefaction  on  the  capillary  throughput,  and 
capillary  end  effects. 

2  Scope  of  Simulations 

In  the  present  work,  two  reservoirs  filled  with  a 
single-species  gas  and  joined  by  a  two-dimensional 
capillary  of  a  finite  length  L  and  height  h  are  con¬ 
sidered  in  the  computations.  The  two-dimensional 
configuration  corresponds  to  the  case  of  a  capillary 
of  infinite  width  w.  A  schematic  of  the  geometric 
configuration  and  notations  used  throughout  this 
paper  are  shown  in  Fig.  1.  The  subscripts  I  and 
II  refer  to  the  quantities  attributed  to  the  left  and 
right  container,  respectively  (far  enough  from  the 
channel  inlet  and  outlet).  The  channel  walls  have 
some  temperature  distribution  varying  from  Tj  to 
Tjj.  If  the  reservoirs  are  at  the  same  temperature 
(Tj  =  Tjj)  and  the  pressure  pj  is  larger  than  pjj 
there  wTill  be  a  mass  flow  from  the  left  container  to 
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the  right  one.  Conversely,  if  the  pressures  are  equal 
but  a  temperature  drop  is  maintained  between  the 
reservoirs  then  the  gas  flow  develops  from  the  cold 
container  to  the  hot  one. 

The  summary  of  cases  that  have  been  considered 
for  two-dimensional  channel  flows  is  given  in  Table 
1.  The  gas  is  molecular  nitrogen  in  all  simulations. 
The  first  set  of  computations  is  for  a  pressure-driven 
flow  in  finite  length  channels  for  a  range  of  Knud¬ 
sen  numbers  from  near-continuum  to  free  molecular 
(cases  I-a  to  I-e).  Hereafter,  the  Knudsen  number 
is  based  on  the  mean  free  path  in  the  right  reservoir 
and  the  channel  height.  Note  also  that  the  change 
in  Knudsen  number  is  due  a  change  in  the  mean 
free  path  and  not  the  channel  geometry.  A  second 
set  includes  cases  Il-a  to  Il-e  that  correspond  to 
a  flow  with  both  pressure  and  temperature  gradi¬ 
ents  in  a  channel  with  a  length-to-height  ratio  of  5. 
For  each  of  these  five  cases  we  consider  two  dif¬ 
ferent  wall  temperature  distributions:  linear  and 
step-wise.  Finally,  the  influence  of  different  wall 
temperature  distributions  is  examined  in  a  channel 
that  connects  reservoirs  with  a  temperature  ratio 
Tj/Tjj  =  2  and  zero  pressure  drop  (cases  Ill-a  to 
III-c). 

2.1  DSMC  Method 

The  DSMC  method  [6]  has  been  applied  in  this 
work  to  obtain  numerical  solutions  for  short  mi- 
crochannel  flows.  The  DSMC  method  is  a  numer¬ 
ical  approach  to  model  gas  flows  in  a  wide  range 
of  Knudsen  numbers  from  the  near  continuum  to 
free  molecular  flow  regimes.  The  DSMC-based  soft¬ 
ware  SMILE  [7]  is  used  for  all  DSMC  computations. 
The  code  uses  the  majorant  frequency  scheme  [8]  of 
the  DSMC  method  for  modeling  of  collisional  pro¬ 
cess.  The  intermolecular  potential  is  assumed  to 
be  the  variable  soft  sphere  model  [9].  The  Larsen- 
Borgnakke  model[10]  with  temperature-dependent 
Zr  and  Zv  and  discrete  rotational  and  vibrational 
energies  is  used  for  the  energy  exchange  between 
the  translational  and  internal  modes.  At  the  gas- 
wall  interface  the  the  Maxwell  model  is  used  with 
the  fully  diffuse  accommodation  of  tangential  mo¬ 
mentum  and  energy. 

2.2  BGK  Model  Kinetic  Equation 
Solution 

The  computational  cost  of  DSMC  simulations  of 
low-speed  flows  in  long  channels  becomes  pro¬ 
hibitive  due  to  low  signal-to-noise  ratio  [11].  This 
is  why  for  modeling  of  such  flows  it  is  reasonable  to 
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consider  the  solution  of  the  steady-state  Boltzmann 
equation  with  the  collision  operator  simplified  using 
BGK  model[12], 

t,s  +  vf  =  (1) 

where  f(x,y,U,V,W)  is  the  distribution  function, 
x  and  y  are  Cartesian  coordinates,  U,  V,  W  are  ve¬ 
locity  components,  v  is  the  collision  frequency  and 
/o  is  the  local  equilibrium  Maxwell-Boltzmann  dis¬ 
tribution  function  which  depends  on  the  moments 
of  the  distribution  function  /.  This  is  a  non-linear 
integrodifferential  equation  that  uses  a  simplified 
form  of  the  collision  integral  compared  to  the  Boltz¬ 
mann  equation,  and  is  applicable  to  model  gas  flows 
with  arbitrary  Knudsen  numbers  and  degree  of  flow 
nonequilibrium. 

The  BGK  model  kinetic  equation  can  be  solved 
using  a  finite-difference  method  and  standard 
discrete-ordinate  method  procedures  [13]  for  inter¬ 
polation  of  the  distribution  function  in  velocity 
space.  The  velocity  components  ( U ,  V)  are  trans¬ 
formed  into  polar  coordinates  (S,<j>).  The  Gauss- 
Hermite  half  range  quadrature  [14]  of  order  16  and 
three-eights  rule  are  used  for  integration  over  the 
speed  and  velocity  angle,  respectively. 


3  Results  and  Discussion 

3.1  Pressure-Driven  Flow  in  Short 
Channels 

First  consider  the  pressure-driven  flow  in  a  short 
channel  with  a  length- to-height  ratio  of  5.  The 
DSMC  solutions  have  been  obtained  for  different 
Knudsen  numbers  from  the  near  continuum  to  free 
molecular  regime  (cases  I-a  to  I-e  in  Table  1).  The 
ratio  of  the  upstream  to  downstream  reservoir  pres¬ 
sures  in  these  cases  is  1.5. 

The  end  effects  in  a  short  channel  result  in  the 
effective  elongation  of  the  capillary  which  is  clearly 
observed  in  Fig.  2  where  the  pressure  profiles  along 
X-axis  are  plotted  for  different  Knudsen  numbers. 
The  pressure  drop  extends  significantly  beyond  the 
channel  inlet  and  outlet  for  all  Knudsen  number 
cases.  The  deviation  of  pressure  from  a  linear  pro¬ 
file  is  shown  in  Fig.  3.  It  is  seen  that  such  a  devia¬ 
tion  increases  with  Knudsen  number  and  reaches  a 
maximum  of  about  6%  for  Kn//=50. 

The  calculated  mass  flow  rates  for  Cases  I-a  to 
I-e  are  given  in  Table  4.  The  flow  coefficients  Qp 
and  Qt  are  used  here  that  are  usually  employed  for 
characterization  of  the  Poiseuille  and  thermal  creep 
flows  for  the  conditions  of  small  temperature  and 


pressure  gradients.  The  flow  coefficients  are  defined 
as  follows: 


Qp 


Qt 


M 

M 

(l^WT0)PoWh^^ 


(2) 

(3) 


where  M  is  the  mass  flow  rate  (kg/s),  po,  Po,  and  To 
are  the  average  density,  pressure,  and  temperature, 
respectively,  and  w,  h,  L  are  channel  width,  height 
and  length,  respectively. 

The  calculated  flow  coefficient  Qp  is  plotted  in 
Fig.  4  as  a  function  of  Knudsen  number.  There  are 
also  shown  the  flow  coefficient  Qp  for  a  Poiseuille 
flow  between  parallel  plates[15]  for  small  pressure 
differences  and  the  flow  coefficient  Qp  for  a  short 
channel  with  L/li=5  in  the  free  molecule  regime 
[16].  As  expected,  the  flow  coefficient  is  larger  for 
the  longer  channel  due  to  the  smaller  influence  of 
end  effects  on  the  mass  flow.  Note  also  that  the 
computed  free  molecular  values  agree  well  with  [16]. 
The  difference  between  the  computed  flow  coeffi¬ 
cients  and  the  one  calculated  in  [15]  becomes  larger 
at  large  Knudsen  numbers.  The  minimum  in  the 
flow  coefficient  observed  at  a  Knudsen  number  of 
about  3  illustrate  the  well  known  Knudsen  mini¬ 
mum  effect. 

The  appearance  of  a  conductance  minimum  ob¬ 
served  by  Knudsen  has  been  qualitatively  explained 
by  Pollard  and  Present  [17].  They  suggested  that  at 
low  pressure  conditions,  when  the  mean  free  path  of 
the  molecules  is  much  greater  than  the  height  and 
the  length  of  the  channel,  there  are  molecules  that 
have  near  axial  velocity  components  on  entry  to  the 
channel  or  after  a  collision  with  the  surface.  These 
molecules  tend  to  make  a  much  larger  contribution 
to  the  mass  flow  than  other  molecules.  When  pres¬ 
sure  increases,  the  fraction  of  such  molecules  de¬ 
creases  since  their  paths  are  disrupted  by  molec¬ 
ular  collisions.  On  the  other  hand,  the  molecu¬ 
lar  collisions  generally  tend  to  create  a  drift  ve¬ 
locity  that  increases  mass  flow.  Since  the  develop¬ 
ment  of  a  drift  velocity  depends  on  the  height-based 
Knudsen  number,  and  the  number  of  molecules  that 
go  straight  through  the  channel  depend  on  a  the 
length-based  Knudsen  number,  the  impact  of  the 
decreasing  molecular  mean  free  path  initially  out¬ 
weighs  the  drift  velocity  effect.  As  a  result,  the  con¬ 
ductance  initially  decreases  with  Knudsen  number, 
and  therefore  creates  a  minimum. 

To  quantitatively  examine  the  Knudsen  minimum 
effect,  DSMC  computations  has  been  performed  in 
a  channel  of  a  length-to-height  ratio  of  15  using 
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five  chemically  identical  species.  Species  1  and  4 
are  assigned  to  molecules  that  come  from  the  up¬ 
stream  and  downstream  boundaries,  respectively. 
When  molecules  of  species  1  and  4  collide  with  the 
channel  surface,  they  are  transformed  to  species  2 
and  5,  respectively.  Additionally,  to  study  the  ef¬ 
fect  of  molecular  collisions,  species  3  is  created  from 
species  1  and  2  when  they  experience  collisions  with 
any  other  species  inside  the  channel.  Species  4  that 
collides  with  other  molecules  inside  the  channel,  is 
transformed  to  species  5. 

The  distribution  of  mole  fractions  of  these  species 
along  the  centerline  for  Kn//=75  is  shown  in  Fig.  5. 
Here,  the  channel  starts  at  X=0.005  m,  and  ends  at 
X=0.020  m.  Note  first  that  although  the  popula¬ 
tion  of  molecules  with  no  surface  or  molecular  colli¬ 
sions  (Species  1)  decreases  significantly  through  the 
channel,  but  still  remains  at  about  1%  at  the  exit. 
Species  1  is  gradually  replaced  by  species  2  and 
3,  with  the  former  gradually  decreasing  along  the 
channel,  and  latter  reaching  its  maximum  near  the 
center.  Note  also  that  near  the  exit,  the  population 
of  molecules  that  come  from  the  upstream  bound¬ 
ary  decreases  to  about  10%  due  to  the  backscatter- 
ing.  The  remaining  90%  of  molecules  come  from 
the  downstream  reservoir.  The  fraction  of  these 
molecules  that  travel  upstream  through  the  channel 
and  reach  the  entrance  is  about  0.3%  (Species  4). 
Consider  now  the  contribution  of  different  species 
into  the  total  mass  flow  given  in  Table  2  along  with 
their  relative  concentration.  The  important  fact  is 
that  although  the  concentration  of  species  1  is  about 
three  times  smaller  than  species  2,  its  contribution 
to  the  mass  flow  is  only  about  two  times  smaller. 
Another  important  fact  is  that  species  1  effectively 
contributes  about  40%  to  the  total  mass  flow,  since 
the  the  contributions  of  species  3,  4,  and  5  cancel 
out.  That  means  that  if  species  1  had  an  X-velocity 
similar  to  those  of  species  2  and  3,  the  total  mass 
flow  would  be  almost  12%  lower. 

The  results  for  a  smaller  Kn//=1.5  are  presented 
in  Fig  6.  The  most  important  conclusion  here  is 
that  no  molecules  cross  the  channel  without  colli¬ 
sions  with  other  molecules.  The  mole  fraction  of 
species  1  and  2  at  the  channel  exit  and  species  4 
at  the  channel  entrance  are  essentially  zero.  It  is 
interesting  to  note  that  the  fraction  of  upstream 
molecules  reaching  the  exit  is  about  10%,  similar  to 
the  larger  Knudsen  number  case.  Comparison  of  the 
species  mass  flows  is  shown  in  Table  3.  Since  the 
mass  flow  of  species  4  increases  50  times  as  com¬ 
pared  to  Kn//=75,  same  as  the  density  does  for 
these  two  cases,  it  is  clear  that  there  is  no  yet  any 
drift  velocity  formed  for  Kn//=1.5.  The  absence 
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of  a  high-velocity  contribution  from  species  1  ex¬ 
plains  therefore  the  Knudsen  minimum,  with  the 
mass  flow  per  pressure  drop  decreasing  when  the 
Knudsen  number  is  reduced  from  75  to  1.5.  The 
computations  also  showed  that  a  drift  velocity  be¬ 
gins  to  form  at  a  Knudsen  number  of  about  0.75. 

The  existence  of  the  Knudsen  minimum  is  illus¬ 
trated  in  Fig.  7  where  the  mass  flow  per  pressure 
drop  is  shown.  The  reduced  mass  flow  is  about 
18%  higher  at  Kn//=75  than  at  Kn//=1.5.  About 
two  thirds  of  this  is  due  to  molecules  that  cross  the 
whole  channel  without  gas  and  surface  collisions, 
and  the  remaining  part  is  due  to  molecules  that  af¬ 
ter  a  collision  with  the  surface  acquire  low  radial 
and  high  axial  velocity  components.  The  compu¬ 
tations  also  showed  that  the  Knudsen  minimum  is 
observed  not  only  for  a  relatively  long  channel  of 
L/h=15,  but  also  for  a  short  channel  with  L/h=5. 
For  the  latter  one  the  minimum  amounts  to  only 
about  1.5%. 

3.2  Pressure  and  Temperature 
Driven  Flow  in  a  Short  Channel 

In  the  above  section,  the  flow  of  a  gas  driven  by 
density  gradients  was  considered  in  the  absence  of 
a  temperature  difference  between  reservoirs.  The 
situation  is  often  reversed  in  a  Knudsen  compres¬ 
sor  nanomembrane  flow,  when  the  flow  is  primar¬ 
ily  driven  by  temperature  gradients.  In  order  to 
study  separately  the  temperature  effect,  a  temper¬ 
ature  drop  is  assumed  across  the  channel,  with  the 
density  in  the  reservoirs  being  equal.  Such  flows 
through  finite  length  channels  in  the  presence  of 
both  temperature  and  pressure  differences  are  ex¬ 
amined  in  this  section.  The  presented  results  corre¬ 
spond  to  Cases  Il-a  to  ITe  in  Table  1.  In  all  of  the 
considered  cases,  the  pressure  ratio  Pi/Pn  =  1.5 
and  temperature  ratio  is  Tj/Tjj= 1.5.  Since  the 
pressure  ratio  in  these  cases  exceeds  the  thermo- 
molecular  pressure  ratio,  there  is  a  non-zero  mass 
flow  from  the  left  reservoir  to  the  right  one.  For 
each  set  of  reservoir  conditions  two  different  cases 
of  temperature  distribution  are  simulated,  (1)  a  lin¬ 
ear  distribution  from  Tj  to  T/j,  and  (2)  a  step-wise 
distribution  Tw[x )  =  Tj  for  x  <  L/ 2;  Tw{x )  =  Tn 
for  x  >  i/2. 

The  calculated  X-component  of  the  velocity  field 
is  shown  in  Fig.  8  for  the  Cases  Il-d  and  ITe  that 
correspond  to  the  Knudsen  number  of  0.2  and  0.05, 
respectively.  The  most  important  conclusion  that 
may  be  drawn  from  these  figures  is  that  the  temper¬ 
ature  distribution  does  not  impact  the  flow  outside 
of  the  channels.  This  is  true  for  both  Knudsen  num- 
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bers  considered.  Moreover,  the  velocity  fields  for 
the  linear  and  step-wise  temperature  distributions 
are  close  near  the  surface,  but  differ  significantly  at 
the  centerline. 

Comparison  of  density  profiles  along  the  center- 
line  for  the  above  cases  is  given  in  Fig.  9.  General 
behavior  of  density  profiles  is  similar  for  different 
Knudsen  numbers  and  temperature  variation,  with 
density  being  lower  in  the  high-temperature  part, 
and  higher  in  the  low-temperature  part.  This  is 
related  to  a  relatively  lower  thermal  velocities  of 
the  latter  and,  as  the  result,  their  longer  residence 
time  and  larger  densities.  Similarly,  the  population 
of  molecules  in  the  high  temperature  region  is  af¬ 
fected  by  short  residence  times  in  that  region.  This 
effect  is  obviously  more  pronounced  in  the  step  dis¬ 
tribution  cases.  The  molecular  collisions  increase 
the  deviation  of  density  from  its  average  value. 

The  X-component  of  velocity  profile  along  the 
centerline  is  plotted  in  Fig.  10.  The  pressure- 
driven  Poiseuillc  and  temperature-driven  thermal 
creep  flows  represent  two  competing  flow  mecha¬ 
nisms  in  these  cases.  Due  to  the  thermal  creep 
flow  from  the  cold  container  to  the  hot  one,  the 
X-component  of  velocity  is  lower  in  Cases  Il-e  (lin) 
and  Il-e  (step)  than  in  Case  I-e.  In  the  case  of 
step-wise  wall  temperature  distribution,  the  result¬ 
ing  velocity  profile  in  the  channel  is  non-monotonic, 
and  has  a  strong  gradient  in  the  step  region  due  to 
thermal  creep. 

3.3  Influence  of  Wall  Temperature 
Distribution 

The  influence  of  wall  temperature  distribution  on 
the  thermal  transpiration  flow  in  short  microchan- 
nel  was  considered  for  a  channel  with  lengtlr-to- 
height  ratio  of  5  at  Knn= 0.2.  To  study  this  effect, 
the  pressure  was  assumed  equal  in  the  two  reser¬ 
voirs,  and  the  temperature  drop  was  set  to  300  K. 
The  three  cases  of  temperature  distribution  are  (i) 
the  linear  temperature  variation  along  the  channel 
wall  (case  Ill-a),  (ii)  a  step-wise  temperature  vari¬ 
ation  (case  Ill-b)  and  (iii)  a  non-monotonic  distri¬ 
bution  typical  for  a  radiatively  heated  nanomem¬ 
brane  Knudsen  compressor  [18]  (case  III-c).  A  non¬ 
monotonic  wall  temperature  profile  has  a  maximum 
of  750  K  at  one-fourth  of  the  channel  length. 

The  change  in  the  wall  temperature  distribution 
significantly  impacts  the  flow  velocity  (see  Fig.  11). 
Note  that  only  a  portion  of  computational  domain  is 
shown  to  illustrate  the  details  of  the  flow  inside  the 
channel.  This  impact  is  noticeably  stronger  that  the 
one  shown  in  the  previous  section,  and  is  observed 
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both  at  the  centerline  and  near  the  surface. 

The  influence  of  wall  temperature  distributions 
inside  the  channel  on  the  centerline  temperature 
profile  is  illustrated  in  Fig.  12.  The  gas  temperature 
generally  follows  the  surface  temperature  inside  the 
channel,  with  the  end  effects  being  visible  even  at  a 
distance  from  the  inlet/outlet  equal  to  the  channel 
length.  The  end  effects  are  somewhat  smaller  for 
the  stepwise  distribution,  and  largest  for  the  non¬ 
monotonic  profile. 

Figure  13  shows  the  calculated  pressure  profiles 
along  the  centerline  for  the  above  temperature  dis¬ 
tributions.  The  qualitative  behavior  is  similar  for 
the  linear  and  stepwise  distribution,  with  pressure 
larger  in  the  hot  part,  and  smaller  in  the  cold  part. 
The  non-monotonic  case  significantly  differs  from 
the  first  two,  since  there  is  no  significant  pressure 
minimum  in  the  cold  part.  Also,  there  is  a  pressure 
increase  inside  the  channel  of  about  10  percent  for 
the  monotonic  profile  due  to  the  reverse  transpira¬ 
tion. 

The  calculated  mass  flow  rate  and  flow  coeffi¬ 
cient  Qt  for  the  three  cases  are  listed  in  Table  4. 
The  mass  flow  rate  in  case  of  stepwise  temperature 
variation  is  only  slightly  lower  (by  «  3%)  than  for 
the  linear  temperature  variation.  However,  the  re¬ 
verse  transpiration  in  the  non-monotonic  tempera¬ 
ture  distribution  case  results  in  significant  degrada¬ 
tion  of  mass  throughput  (by  «  18%). 

3.4  BGK  Model  Solution  for  a  Long 
MicroChannel 

Results  of  the  validation  study  of  a  long  channel 
flow  modeling  using  a  BGK  model  kinetic  equation 
are  presented  in  this  section.  The  flow  of  nitrogen 
in  a  channel  with  a  length-to-height  ratio  of  30  was 
calculated  for  a  pressure  ratio  of  1.5  and  a  constant 
temperature  of  the  wall  of  300  K.  The  tempera¬ 
ture  in  the  upstream  and  downstream  reservoirs  is 
equal  to  300  K.  The  Knudsen  number  based  on  the 
channel  height  and  the  mean  free  path  in  the  down¬ 
stream  reservoir  is  0.05. 

The  X-component  of  velocity  field  is  shown  in 
Fig.  14  obtained  using  three  different  uniform  spa¬ 
tial  grids  in  the  BGK  solver.  The  grid  dimensions 
are:  (i)  101x51,  (ii)  51x26  and  (iii)  26x14.  The  an¬ 
alytical  approximate  solution[19]  of  Navier-Stokes 
equations  for  a  long  channel  flow  with  slip  bound¬ 
ary  conditions  is  also  plotted  in  Fig.  14.  The  discrete 
ordinate  solution  of  BGK  model  kinetic  equation 
with  coarser  grids  results  in  a  lower  velocity  in  the 
whole  computational  domain.  Note  that  a  similar 
dependence  of  results  on  the  number  of  simulated 
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particles  and  cells  in  the  DSMC  method  was  ob¬ 
served  [20].  The  solution  obtained  for  a  101x51  grid 
is  in  excellent  agreement  with  the  Navier-Stokes  so¬ 
lution. 

4  Conclusions 

The  direct  simulation  Monte  Carlo  method  and  the 
discrete  ordinate  method  for  the  solution  of  the 
BGK  model  were  used  to  numerically  study  rarefied 
gas  flows  through  two-dimensional  finite  length  cap¬ 
illaries.  A  wide  range  of  Knudsen  numbers  was  ex¬ 
amined,  from  0.05  to  75,  in  the  presence  of  den¬ 
sity  and/or  temperature  gradients.  The  lengtli-to- 
height  ratios  varied  from  5  to  30.  Three  differ¬ 
ent  wall  temperature  distributions  were  considered, 
namely,  linear,  step-wise,  and  a  non-monotonic  pro¬ 
file  typical  for  a  radiantly  heated  Knudsen  compres¬ 
sor  membrane. 

For  the  pressure-driven  flow  in  a  short  channel 
with  a  length-to-height  ratios  of  5  and  15,  the  end 
effects  result  in  the  effective  elongation  of  the  cap¬ 
illary.  The  effective  elongation  increases  with  the 
Knudsen  number.  The  flow  coefficient  is  larger 
for  the  longer  channel  due  to  smaller  influence  of 
end  effects  on  the  mass  flow.  The  minimum  in  the 
flow  coefficient,  so  called  Knudsen  minimum,  is  ob¬ 
served  at  a  Knudsen  number  of  about  3  for  a  chan¬ 
nel  with  L/h=15.  The  mechanism  for  the  Knudsen 
minimum  in  flow  conductance  is  quantitatively  ex¬ 
plained  by  using  five  chemically  identical  species  in 
DSMC  simulations.  The  computations  showed  that 
the  Knudsen  minimum  exists  not  only  for  a  rela¬ 
tively  long  channel  of  L/h=15,  but  also  for  a  short 
channel  with  L/h=5. 

The  influence  of  wall  temperature  distribution  on 
the  thermal  transpiration  flow  in  a  short  channel 
at  Knu= 0.2  was  considered.  The  qualitative  flow 
behavior  is  similar  for  the  linear  and  stepwise  dis¬ 
tribution,  with  pressure  larger  in  the  hot  part,  and 
smaller  in  the  cold  part.  The  non-monotonic  wall 
temperature  distribution  case  significantly  differs 
from  the  first  two,  since  there  is  no  significant  pres¬ 
sure  minimum  in  the  cold  part. There  is  a  pressure 
increase  inside  the  channel  of  about  10  percent  for 
the  monotonic  profile  due  to  the  reverse  transpira¬ 
tion.  The  mass  flow  rate  in  case  of  stepwise  temper¬ 
ature  variation  is  only  slightly  lower  (by  3%)  than 
for  the  linear  temperature  variation.  However,  the 
reverse  transpiration  caused  by  the  non-monotonic 
wall  temperature  distribution  results  in  significant 
degradation  of  mass  flow,  by  about  18%. 

The  application  of  the  numerical  techniques  dis¬ 
cussed  in  this  paper  to  detailed  investigations  and 
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optimizations  of  microscale  thermal  creep  flows 
will  be  extremely  valuable.  Temperature  gradient 
driven  flows  are  very  sensitive  to  thermal  boundary 
conditions.  For  practical  microdevices  it  is  gener¬ 
ally  difficult  to  attain  idealized  thermal  boundary 
requirements.  Thus,  numerical  methods  for  reliably 
estimating  the  magnitude  of  the  consequences  in¬ 
duced  by  non-ideal  conditions  are  important.  The 
work  reported  here  is  a  first  look  at  describing  prac¬ 
tical  thermal  gradient  driven  or  thermal  creep  flows. 

As  an  illustration  of  the  usefulness  of  this  ap¬ 
proach  consider  the  recent  development  by  Young 
[5]  of  a  prototype  Knudsen  Compressor.  Radiant 
energy  was  used  to  generate  temperature  gradients 
in  the  carbon  doped  aerogel  membranes  that  are 
the  active  elements  in  the  compressor’s  pumping 
stages.  Due  to  finite  photon  mean  free  paths  in 
the  carbon  doped  aerogel  membranes,  the  wall  tem¬ 
peratures  in  the  membrane’s  flow  passages  had  a 
maximum  near  to,  but  not  at,  the  exit  surface  of 
the  membrane.  This  was  the  motivation  for  the 
calculations  reported  in  Sec.  3.3,  with  the  non¬ 
monotonic  wall  temperature  variation  of  case  III-c 
constructed  to  mimic  the  essentials  of  a  membrane 
in  Young’s  experiments.  While  the  available  results 
are  at  present  both  scarce  and  preliminary,  they 
do  validate  in  this  case  what  was  anticipated,  and 
provide  useful  quantitative  estimates  of  the  conse¬ 
quences  of  the  particular  wall  temperature  maxi¬ 
mum  chosen  for  the  calculations. 
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Table  1:  Summary  of  cases  considered 


Flow  Conditions  Designation 


Pressure-Driven  Flow  in  Short  Microchannels 

L/h=5,  Na,  77 

=  Tn  =  300  K,  ^=1.5 

Kn  n=50 

case  I- a 

Knu= 3 

case  I-b 

Knn= 1 

case  I-c 

Knii—  0.2 

case  I-d 

Ann- =0.05 

case  I-e 

Pressure  and  Temperature-Driven  Flow 

L/h=5,  N2,  Tv 

=600  K,  T//=300  K,  ^=1.5 

Knu=50 

case  Il-a 

Knii— 3 

case  II-b 

Knu=l 

case  II- c 

Knn=  0.2 

case  II-  d 

Knn= 0.05 

case  II-  e 

Effect  of  Wall  Temperature  Distribution 

L/h=5,  Kn//= 

0.2,  N2,  17=600  K,  77/ =300  K,  ^=1.0 

linear 

case  I I I- a 

step 

case  Ill-b 

non-monotonic 

case  1 1  I-c 
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Table  2:  Average  exit  mole  fraction  and  mass  flow  Table  5:  Calculated  mass  flow  rate  and  flow  coeffi- 
of  different  species  for  Kn//=75.  dent  for  temperature-driven  flow  in  short  channel. 


Species 

Fraction,  % 

m,  kg/s 

1 

1.1 

3.56-10-9 

2 

3.5 

8.14-10-9 

3 

6.2 

1.438-10'8 

4 

46.2 

-9.58-10”8 

5 

43.1 

7.84-10"8 

Total 

100 

4.44-10-9 

Case 

to  (kg/s) 

Qt 

Ill-a 

1.580E-6 

0.1189 

Ill-b 

1.532E-6 

0.1153 

III-c 

1.292E-6 

0.0972 

Table  3:  Average  exit  mole  fraction  and  mass 
of  different  species  for  Kn//=1.5. 

Species 

Fraction,  % 

to,  kg/s 

1 

0 

0 

2 

0 

0 

3 

8.9 

7.22-10~7 

4 

35.3 

-3.88-10"6 

5 

55.1 

3.532-10-6 

Total 

100 

3.74-10"7 

L 


Figure  1:  Schematic  of  the  problem  and  notations 


Table  4:  Calculated  mass  flow  rate  and  flow  coeffi¬ 
cients  for  pressure-driven  flow  in  short  channel. 


Case 

771 

(kg/s) 

Qp 

^  x  10Y 

(  kg/s-  Pa) 

I- a 

0.261xl0"7 

1.00 

4.76 

I-b 

0.429xl0"6 

0.989 

4.69 

I-c 

0.134x  10-5 

1.03 

4.88 

I-d 

0.936Exl0~5 

1.44 

6.82 

I-e 

0.746x  10-4 

2.87 

13.6 

Il-a  (lin) 

0.119xl0"7 

0.533 

2.17 

Il-a  (step) 

0.119xl0-7 

0.533 

2.17 

Il-b  (lin) 
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Figure  2:  Normalized  pressure  along  X-axis. 
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Figure  3:  Deviation  from  the  linear  pressure  distri¬ 
bution  along  X-axis. 
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Mole  Fraction 


Figure  5:  Distribution  of  the  mole  fraction  of  differ¬ 
ent  species  along  the  flow  centerline  for  Kn//=75 


Figure  4:  Flow  coefficient  Qp. 


Mole  Fraction 


Figure  6:  Distribution  of  the  mole  fraction  of  differ¬ 
ent  species  along  the  flow  centerline  for  Kn//=1.5 
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„.  _  ,,  n  ,  .  Figure  9:  Normalized  density  profile  along  the  chan- 

ligure  7:  Mass  now  rate  per  pressure  drop.  ,  .  ,  ,  ,  .  ,  T.  . 

nel  lor  different  temperature  variations  and  Knud- 

sen  numbers. 
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Figure  8:  X-component  of  velocity  field  for  Cases 
Il-d  and  Il-e. 


Figure  10:  X-component  of  velocity  profile  along 
the  channel  for  different  temperature  distributions 
and  Knudsen  numbers. 
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Figure  11:  X-component  of  velocity  field, 

Kn//=0.2,  L/h=5  for  different  temperature  varia¬ 
tions. 


Figure  13:  Pressure  distribution  along  the  channel, 
Kn//=0.2,  L/h=5  for  different  temperature  varia¬ 
tions. 


Figure  12:  Temperature  distribution  along  the 
channel,  Kn//=0.2,  L/h=5  for  different  tempera¬ 
ture  variations. 
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Figure  14:  X-component  of  velocity  field, 

Knr/=0.05,  L/li=30. 


